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A kinetically stabilized axially symmetric tandem mirror (KSTM) uses the momen- 
tum flux of low-energy, unconfined particles that sample only the outer end-regions 
of the mirror plugs, where large favorable field-line curvature exists. The window 
of operation is determined for achieving MHD stability with tolerable energy drain 
from the kinetic stabilizer. Then MHD stable systems are analyzed for stability 
of the trapped particle mode. This mode is characterized by the detachment of the 
central-cell plasma from the kinetic stabilizer region without inducing field-line bend- 
ing. Stability of the trapped particle mode is sensitive to the electron connection 
between the stabilizer and the end plug. It is found that the stability condition for 
the trapped particle mode is more constraining than the stability condition for the 
MHD mode, and it is challenging to satisfy the required power constraint. Further- 
more a severe power drain may arise from the necessary connection of low-energy 
electrons in the kinetic stabilizer to the central region. 
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I. INTRODUCTION 

The tandem mirror magnetic-fusion confinement system is a nearly-cylindrical solenoid 
terminated by a set of plug cells. These plugs consist of simple axisymmetric mirror fields. 
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The earliest stability theories [T] predicted that symmetric mirror machines, including tan- 
dem mirrors would be MHD unstable; experiments in the 1980s on both TARA [2] and 
PHAEDRUS [5 [5 tandem mirror facilities confirmed that MHD instability occurs. Re- 
cently, a new stability innovation, the kinetic stabilizer, has been proposed to stabilize 
an axisymmetric tandem mirror. However, the mechanism by which MHD stabilization is 
achieved in the kinetic stabilizer may make the system susceptible to a rapidly growing 
trapped particle instability |5]. Here we investigate this possibility. 

The kinetic stabilizer is a design proposed by R. F. Post p], inspired by the work of 
D. Ryutov and experimental evidence from the Gas Dynamic Trap (GDT), a single mir- 
ror experiment at Novosibirsk [7j p]. Ryutov established theoretically that an otherwise 
MHD-unstable plasma confined between symmetric mirrors can be stabilized if there exists 
sufficient momentum flux from the effluent plasma on the expanding positive curvature field- 
lines outside the mirrors. The momentum flux generalizes the role of the pressure tensor that 
appears in standard MHD theory. This technique of stabilization was experimentally con- 
firmed using the axisymmetric GDT [9J. The GDT operates in a high-coUisionality regime in 
order to keep the loss-cone full. Ryutov noted that the effect of plasma momentum flux flow- 
ing out of the ends of the machine in the positive-curvature expanding-field region outside 
the mirrors is sufficiently strong to overcome the destabilizing curvature contribution from 
the central part of the plasma. This stabilization has been experimentally demonstrated 
[lUj . It has been shown that MHD stability persists to a moderately high-beta value /3 ~ 0.3 
arising just from the exiting momentum flux. 

In the GDT the loss-cone is filled by relatively strong collisions in the central plasma 
region; with sufficient effluent flow, MHD stability is established. However, in the typi- 
cal tandem mirror the ends of the machine are designed to conflne a long mean-free-path 
plasma which only produces a weak effluent. Thus there is a negligible momentum flux in 
the expander region. As a result, the stabilization mechanism designed for the GDT, i.e. 
stabilization from the momentum flux arising from the effluent plasma, does not simply arise 
in a tandem mirror of conventional design. The kinetic-stabilizer concept proposes to solve 
the problem of low effluence by using external ion-beams injected axially into the machine. 
These kinetic-stabilizer ion-beams are injected at small pitch-angles to the magnetic fleld 
so that the beam can propagate towards the higher magnetic field and then refiect before 
reaching the principal confinement region of the mirror plasma. The ion beam then transits 



3 



out of the machine. While entering and exiting the ends of the machine, the beam forms 
an unconfined plasma with an enhanced momentum flux in a localized region of favorable 
curvature. 




FIG. 1: Axisymmetric magnetic flux surfaces in the KSTM. 
The proposed KSTM reactor is a simple axisymmetric tandem mirror. The plasma in 
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the plug regions possesses higher density and energy than the central-cell plasma, producing 
ambipolar traps in the plugs. The three-dimensional structure of the magnetic flux tube 
can be viewed in figure [TJ Figure [2] provides a schematic diagram of the axial magnetic field 
of a KSTM; figure [3] provides a conceptual sketch of the density profile associated with the 
kinetic stabilizer. At the outer wall, on the left-side of the figure figure [3} the magnetic 
field is at its lowest value B = 5w The density of the beam initially rises in the expander 
region as we consider regions inward from the wall. The beam can either be focused on a 
target in the region around B = B]^s or unfocused. In the focused case, the density will 
peak steeply around B = B]^s and then plummet to a low value just beyond B = By^s- The 
focused kinetic stabilizer ion beam is designed to reflect around the region of i?ks which 
bounds the region of good curvature that exists between the target and the wall. In the 
unfocused case the density will continue to rise until most of the kinetic stabilizer beam 
has been reflected by the rising magnetic field. The remnant beam that reaches this region 
possess a nearly spatially-isotropic distribution, which produces a plateau in density. In the 
plateau region, the density of the unfocused beam is appreciably higher than for a focused 
beam. The electron temperature between the wall and the plateau region will likely be 
lower than the ion energy in the injected kinetic stabilizer. However, the temperature of 
the electrons escaping from the central cell will be from 10-100 keV, and thus the ambipolar 
potential associated with the effluent flux will be extremely high. We expect a transition 
region around a stand-off position, B = B^t, where electron temperature transitions from a 
characteristic temperature of the KS region, to the electron temperature of the core plasma. 
In this region the plasma density will rise. 
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FIG. 2: Axial magnetic field B, in the right half of the proposed KSTM reactor. The field-line 
curvature is negligible to the left of the dashed line. The favorable curvature of the expander region 
lies to the right of the dashed line. A solid black line marks the center of the mirror plug. 



The density profile and the relation between density the electric potential will be dis- 
cussed further insection |Tl] where we present the details of the model used for the KSTM. 



In section III we derive MHD stability relations for the KSTM. Insection IV the trapped 
particle instability is examined in the context of a KSTM and insection |V] we summarize 
the salient conclusions of this work. 
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FIG. 3: A conceptual sketch of the density m the expander region of the KSTM. The shape of the 
density at the wall the entrance to the region of favorable curvature -Bks, plateau, stand-off 
point Bst, and the maximum of the plug are shown. 

II. KSTM MODEL 

A. Distribution of ions in the kinetic stabilizer beam 

The kinetic stabilizer design consists of an ion beam that is injected from outside the 
machine into the expander region of the tandem mirror where the curvature is positive. The 
beam of ions is aligned nearly parallel to the magnetic field. The small pitch-angle of the 
particles in the beam corresponds to injection of ions with small magnetic-moment-per-unit- 
mass /i compared to Eq/B (where Eq is the mean injected-beam-energy-per-unit-mass). We 
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choose a distribution function for the ions in the kinetic stabihzer beam such that they 
possess a fixed energy equal to Eq. The magnetic moment is centered around a value /ig 
with a narrow range of magnetic moments A/x. The distribution function, F, is chosen to 
be 

F(E.,)^Im.(E-E„)^_^^. (1) 

Here the quantity F'j^^^ is related to the particle flux per unit magnetic flux Fks from a single 
end of the tandem mirror as indicated in Q. The spread in magnetic moments of the kinetic 
stabilizer beam is determined by Afi; we use /ix = Eq/Bt^ where the subscript T refers to 
the target position; the values of AyU and /ix effects the maximum density at the target. 
We neglect the ion electrostatic potential energy under the assumption that the potential 
energy is proportional to Tgks, the electron temperature in the kinetic stabilizer region. The 
electron temperature is assumed to be significantly less than the ion beam energy. 

We will treat two types of kinetic-stabilizer ion-beams. In the case of a focused beam, the 
beam is focused on a target outside of the plug /^t- For the second type of beam, /it = 0; 
this will be referred to as the unfocused case, equivalent to the condition /ix ^ A/i. 

The density and pressure in the outer regions of the KSTM are established by an incoming 
particle flux at the wall. The subscript W refers to the wall position, i.e. z = and r = rw 
In terms of the distribution, the particle flux at the wall vrr^Swrks, is 

vrr^^wFks = vrr^ j dE j dnBy^F{E, jj) , (2) 



Eo/Bw 

We find 

oo 



' ^ ^^^^ ■ 'tan-f^U J^). (4) 



ks 

(/i - /ix)' + A^2 2 ' TT \AfiJ - \Afl 







For the unfocused case, /ix/AyU ^ 1, x ^ 1/2, while for the focused case, /ix/A/i ^ 1, 
X ~ 1- 
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B. Density and pressure in the expander due to the kinetic stabilizer beam 



When the effect of the ambipolar potential on the ions is neglected, the density n{B) is 

E/B 



(5) 





We define the quantity x = ixB/Eq and use the distribution function in ([T]) to express the 
density integral as 



n{B) 



V2T'^B fAi2B\ f dx 



7rVE~o [ Eo J J y^l^ix-B/B^y + iAfiB/Eo)'' 



(6) 
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The inequality 1 ^ B/B^ — 1 is satisfied for all B/B^ < 1 when A^Bt/Eq « 1. The last 
approximation is also accurate when IB/B^ — 1| ^ AfiBT^/Eo, but still small. 

The inequality AfiB/Eo ^ 1, determines a region where the density asymptotes to a 
constant value. We call this region the plateau region, where the density npia is 
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(7) 



Here nw = n^B^^) is the density at the wall. 

For the focused case the density peaks near the target position B = B^. At the target 
the density is given by: 



n{B^ 



1/2 



ny/B^ f Eq 
^/2E^ ^AjiB^ J 25w yAfiBr 
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(8) 



For the unfocused case the density increases monotonically from the wall position to the 
plateau region. We can relate the density at the wall to the density at the entrance to the 
region of positive curvature, which we will call the kinetic stabilizer position, denoted by 
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the subscript ks. Assuming that the axial speed of the beam is small compared to the local 
thermal spread of its speed, i.e. /^iiB^^/Eq <^ 1, then 

The plateau region arises when Eq. By focusing the ion beam to a position B — Bt, 

the peak density in the kinetic stabilizer region is enhanced compared to the unfocused 
case by an approximate factor of {Eq/ AiiBtY^'^. In the plateau region, the density of the 
unfocused case is larger than the focused case by a factor of (E'o/A/xSt)^- 

In magnetized mirror-confined plasmas the pressure is typically anisotropic with different 
values for the pressures p±, and respectively. The mass of ions in the kinetic stabilizer 
beam is taken as mi. For our choice of distribution function, the pressures are 

p{B) = p\i{B)+p^{B) = n{B)m,Eo+p\i{B)/2 , (10) 

oo E/ B 

p\l{B) = 2V2miJdE J d^B^E- ^B F{E,n) , (11) 





oo E/B 



P±iB) = 2m,JdE J -7==^/iE,^,) = n{B)m,Eo - ^ . (12) 







The integral for p\\ (B) is 



1 



„,IB) = E^TOiFLs/^ / ^ 5 , (13) 

^2,„,r;.v/2^(i-|)"%(i-|) , (14) 



X 

B-_ 

1/2 

-> 2m^rL\/2E~n 

where 9{x) is the Heaviside step function. For the focused case where AiiB^/Eq <^ 1, the 



perpendicular pressure peaks near the target, B ^ B^, and is approximately 

p^B) ^ m,E^n{B) , (15) 

while the parallel pressure remains small. Thus the ratio of pressure at the target to that 
at the wall is proportional to the ratio of magnetic fields at these positions 

1/2 



(16) 
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In the plateau region the pressure is isotropic and constant with 



p^{B)=p\\{B) = -V2 



TT 



(17) 



3 

V Eo 

When Ayui?T ^ Eq, which is the definition of the unfocused case, the pressure in the plateau 
region is p±{B) = p\\{B) = \/2tt{A/3){T'^^/\^Eo){Eq/ Afi). This pressure is larger than the 
focused case by a factor of {Eq/ AfiB^)^. 

III. KSTM MHD STABILITY RELATIONS 

During steady-state operation, the power P required to sustain the power drain of a 
single-pass kinetic-stabilizer beam is 

P = 27rm,Eor^^o • (18) 

This expression uses the magnetic fiux ipQ = I-Bq'^o center of the central-cell of the 

tandem mirror. The limiting condition for achieving break-even for reactor engineering is 
that the power drain from the kinetic stabilizer beams should not exceed the alpha particle 
energy production rate of the central-cell by more than some factor i]: 

r)nrlm,EoBoT^. < ^Ip^M . (19) 

where rj = 1 is the natural choice. The temperature in this expression is a sum of electron 
and ion temperatures Tc = Tec+Tjc and the subscript c refers to the central-cell. The Lawson 
criterion (nr) f^^s = 2 x 10^^ should hold for our model KSTM. For MHD stability, we 
require 

J dz (pii +p^) > J dz (pii +px) r^^. (20) 

ks plug + central cell 



To evaluate the MHD stability criterion in (20) we constrain the field line radius r(z) to have 



as large as possible variation within the kinetic stabilizer region that is compatible with the 
paraxial approximation. Thus, we take field-line radius r{z) must satisfy 

(fr(z) , , d f dr(z)\'^ (Td 



dz^ dr \ dz ) 2 



We will refer to (21 ) as the marginal paraxial constraint. For optimal stability, the constant 



dp in the the marginal paraxial constraint should be chosen to maximize the integrand on the 
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left hand side of (20), subject to the vahdity of the paraxial approximation for the field-line 



curvature; in this work we use <7p = 1 while acknowledging that additional study is needed 
to determine an optimal value. 



Solving (21) for r{z) yields 
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(22) 



When dr{z-ks)/dz = 0, the inequality in (22) is an equality. Since the choice of dr{z]^<^)/dz = 
maximizes the favorable MHD response, this field-line slope is used in our model for KSTM 
field-lines. 



A. Stability integral in the plug 

We model the field-line radius in the plug region with the form 

r(z) = r.,„„ (ii^El + cos {'^^) ) . (23) 

where Xpi = rmnp/^mxp- The radii rmxp and r^^ip are the maximum and minimum radii of 
the plug region, and < \z — < Lpi. We define the composite pressure p = {p± +p|j)/2. 
We also define 

Pp\ = 777— / dz {p^ + = — — I dzp . (24) 



2Lpi J 2Lpi 



pi 



The dimensionless plasma parameter in this region is defined (3^^ = 2fiQp^i/ B^^^, where 
-Bmnp is the magnetic field at the minimum field-line radius in the plug region. We define /pi 
to be the stability integral over the plug region for the particular magnetic field assumed in 



(23). Thus we need to evaluate 

d'^r 



/pi = -2 / dzpy— . (25) 



pi 
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We indicate an approximation of Jpi with a prime. We evaluate /pj, as the approximate 
stabihty integral Jpi when the pressure is taken as isotropic in the plug with a pressure Ppj. 
We obtain 



= ^^u(i-pini+7(^) I- (26) 




In this expression we have used the paraxial approximation for the magnetic flux ?"mxp-^mnp ~ 
Bor^. For simphcity we have neglected the small factor (1 — Xpi)^/(2 (1 + Xpi))^ in the 
estimate /pj. We define Opi as a dimensionless, order- unity parameter that is the ratio of the 
exact to the estimated value of the stability integral 

«pi = (27) 

The factor a^i may be less than unity {e.g. due to the inclusion of compressibility effects 
in the plug cell and improved design of the plug, etc.) and such improvement gives some 
flexibility in estimating a window of operation for the KSTM. 

The contribution to MHD stability from the central-cell must also be considered and 
is treated in the same way as the plug. We use a shape of the same form as as the plug 



region described in (23), and central-cell parameters rmxp — ^ '"O; -^pi 2Lc, z — Lc ^ z, 
Xp\ Xc = ro/rj, 

(28) 

The central-cell stability contribution will be used below. 




B. Stability integral in the expander region 

For our KSTM model there is no contribution to the stability integral from the inter- 
mediate region between the outlet of the plug (defined by the maximum of the magnetic 
field) and the beginning of the positive-curvature kinetic stabilizer region at z = z^g. The 
stability contribution of this intermediate expander region would vanish if either (Pr/dz"^ or 
the kinetic pressure were negligible. 

The MHD contribution of the stability drive integral from the expander region is 

/ks = j in + p^) ^ ■ (29) 
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In this expression k is the field-hne curvature, 



1 



/drV 
(d~z 



d'^r 



3/2 dz^ 



(30) 



The paraxial approximation of the curvature is given in the right-most term of (30). 

To estabhsh a base-case for scahng our results, we consider the unfocused kinetic stabilizer 
beam (AyU = 0), and use the paraxial approximation for the curvature. Using the fact that 
magnetic flux is conserved, Br^ = BqTq, the expressions for the pressures from (11) and 



(12), and the solution r{z) from (22), we find 
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where we have used the slope dry^/dz = 0. We approximate the stability integral in the 
expander region as 



d^ 



(31) 
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= 2mi^/o~2E^ BorlT'^^r^Jv ■ 



We define the ratio of the exact to the approximate stability integrals 
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(32) 
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In (31) the use of the paraxial approximation for the curvature allows /ks to be made ar- 
bitrarily large by increasing dp. However, as increases, the paraxial approximation is 
eventually violated. To verify consistency between /ks and the approximate paraxial ap- 
proximation, we choose o"p = 1 and compare results for the non-paraxial expression for the 
field-line curvature with the results of the paraxial approximation. We use the paraxial 



solution for z from (22) and the more accurate form of the curvature k to calculate the 
integral 

zw 1 

rlo'TK r flu 1 

■ (34) 

»'ks/''W 
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For the unfocused case (31) becomes 
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FIG. 4: The quantity a^g = -^ks/-^ks evaluated for an unfocused kinetic stabihzer beam from (35) 
using (Tp = 1. 



Figure |4] shows the comparative accuracy of the paraxial approximation and the more 
accurate non-paraxial curvature for a^g = /ks/-^ks ^ function of x = rks/Tw The maxima 
are ctks = 1-08 for Xks = 0.43 for the paraxial form and Oks = 0.815 for Xks = 0.63 for the 
non-paraxial form. There is some quantitative discrepancy, which indicates that cTp = 1 may 
be the largest value that can be used in the expression for Jks to give reliable results. 

Larger values of Jks can be obtained by focusing the beam so that the target is at the 
kinetic stabilizer point -Bx = -Bks- When the beam is focused on this location, a logarithmic 
divergence arises in Jks, producing arbitrarily strong MHD stabilization. For a beam with 
thermal spread (finite A/i), the magnitude of the stability integral Jks is bounded. To 
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evaluate the stability integral in this case we use an approximate form for the pressure: 

2p 
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(36) 



This is valid if 1 ^ BjB^ — 1, a condition that is always satisfied when the right-hand side 
is negative and /\^Bt/Eq ^ 1. Again using the condition that magnetic flux is conserved, 
so B/B]^s = {fks/^Yi and setting the target position T at the kinetic stabilizer position KS, 
we obtain a stability integral Jfo = /ks for the focused case 



'fo 



ks 9 



(37) 




(38) 



A plot of Oks = I'lo/^Ls ^ function of rks/Tw is shown in figure [s] for the focused case 
where (A/ii?ks/-Eo = 10"^,10~^). As AfiBi^s/ Eq decreases, the optimal MHD response is 
found closer to the wall; at small AfiB]^^/ Eq, the value of this optimal response is given by 
0.2 + 0.71ogio(^o/A/u5ks). For A/iSks/^o = 10"^ this leads to Oks ~ 3. 
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FIG. 5: The quantity a^g = /ks/-^ks focused case, evaluated at two different beam widths 



[(31)]. 



C. Operating regime scaling laws 



For our machine to be stable we require Mks = (/pi + /c) //ks < 1 or 



1 > Mk, 
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where 



^o/^mxp) /c 



(39) 



Tec + Tic. We shall call Mks the MHD marginality parameter. 



The stability condition of (39) uses the approximations for the stability drive integrals in 



the plug given by (26) and in the kinetic stabilizer region given by (31). The MHD drive 



from the central-cell uses the same functional forms for the magnetic field's axial variation 
as the plug region, and employs the appropriate changes in lengths and mirror ratio. The 
total length of the central-cell region is 2Lc. In this stability criterion the central-cell drive 
is integrated over only half of the central-cell. The quantity Oc is the ratio of the actual 
MHD drive to the scaled MHD drive in the central-cell. 
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From the power-balance relation (19) we find 
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Thus the condition that the system be MHD stable with an acceptable energy loss from the 
kinetic stabilizer beam is 
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The plasma beta is f3c = 11X^2^0/ Bq. We use 

_2i Te(100 keV) / 3 



nZHm.-^) = 4.5 ■ 10" 



/3c 



(42) 



Normalizing (42) to a reference parameter of approximately 100 MW fusion power produc- 
tion, we find 

2 
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Bo AeV^o(keV) / Te(100keV) \ / 3 



Here Ac is a dimensionless parameter that characterizes the machine 



(43) 



A. 



1/2 



(44) 



lOapi {1 - ^ + 2.50a^^ {1 - xiy 

The injected ion mass mi (a) in the kinetic stabilizer beam is in atomic units so that m; = 2 
for deuterium. The Lawson number {nT)f^s = 2 x 10^'^ m~'^ is used to characterize the alpha 
particle power production; the terms in each of the parentheses is estimated to be order 
unity and thus Ac — 1 is a characteristic estimate. 
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The scaling law (43) indicates whether a system with a given set of parameters can be 
MHD stable with an acceptable energy loss from the kinetic stabilizer beam. We choose a 
reasonable set of parameters to test the MHD stability: 

fio(T) = 3T, Tc(lOO) = 100 keV, tq = ^ m, 



m / - - ^ 

5 



/3pi = = Oks = ftc = 2Q;pi = 1, 

„2 



= 0.5, xi = 0.1, mi(a) =2, r] = l. (45) 

With this choice of parameters we find Ac = 0.175. We find it necessary to have -Bw < 66 
gauss when Eq = 1 keV or < 660 gauss when Eq = 100 eV. For the parameters 



in (45) the instability drive of the central-cell is 0.4 times the drive from the plug. If 
the drive from the plug could be reduced or eliminated, the drive in the central-cell alone 
would require < T.Q/^^IO^ gauss for = 1 keV (or 0.79P^ T if = 100 eV). Hence 
significant improvement could be achieved if the drive from the plug could be neutralized. 
For example, the plug drive could be neutralized if the particle species in the plug were 
sufficiently energetic |TT] . 

In summary we find that an expander whose magnetic field at the wall could be as large 
as the order of several 100 gauss would enable MHD stability to be achieved by kinetic 
stabilizer beams of several hundred eV, in combination with 3 T central-cell magnetic field 



and ctpi ~ 0.5. The detailed scaling is indicated in (43) and (44) 



D. Beta Limitation and Adiabaticity Limits 

An additional facet of the kinetic stabilizer ion plasma beta should be considered. When /3 
becomes larger than unity, typically the field-lines cannot coUimate the plasma because self- 
consistent effects drive the plasma across field-lines [12]. For the injected kinetic stabilizer 
ions, the local /3ks = 2/iop/-B^ must be less than unity along the imposed field-lines. From 
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(36) we find that in the kinetic stabilizer region 
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(46) 



Expressing FJ^g in terms of /3ks, we find that the MHD stability condition from (39) and (46) 
can be expressed as 
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Here we again use the conservation of flux {rQ/r^){Bo/ B) = (B^^/B). If the inequality has 



to be satisfied at every position, (47) becomes 

1 > Mk« > 



/0"n 



(1 - xiYPcgMx 
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1 + 



4api /3pi (1 - x 



2 \2-\ 



,^2^2 



ttc ^pi /3c (1 - 

Here the expression Mx is the maximum with respect to B of the function 



(48) 
(49) 



Mx = Max 



B 



w 



5 



1/2 



+ V2 1 



"ks ksJ 



(50) 



For the unfocused case where 6t = (also nearly applicable for the focused case 
if eks-Bks/-^w — 1)' ^he maximum of this function occurs at i? = Byj, so Mx = 



2v2a/1 — -Bw/-Bks- For the focused case, when eks-Bks/-^w ^ ^ maximum occurs when 
&ks = &T = 1 - eks/ys, yielding Mx = 1.26(Sw/5ks)(l/V^). 

The best MHD condition is achieved by making Mks as small as possible subject to the 
condition that /3ks always remain less than unity. To fulfill MHD expectations and the 
/3ks < 1 condition limits the allowable range for Mks in the unfocused case (or a focused case 
where Afj,B\^s/Eo > 1) to 



1 > Mks = 0.95^1 - B^/Br^g . 



'Or. 



(51) 
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For the focused case the permitted range of operation for M^s is 

1>M,.>0.4|^^^,. (52) 

Thus the condition that the beta value of the kinetic stabihzer beam be less than unity may 
limit the utility of obtaining enhanced MHD stabilization by focusing the beam at i? = 



because the inequality in (52) could be violated for small eks- 

Another possible limitation to the KSTM is that the upper ratio of i?o/-Bw = ^w/^o ^lust 
satisfy the particle orbit adiabaticity condition. The ion parallel velocity be significantly 
larger than the curvature drift velocity, requiring that 

^/2Eo (T,d > .X , 2 = V ' ^^^^ 

where cXad ~ 1 is a constant associated with the adiabaticity. We find 

O^E(keV) Zap y/mi(a)\ 
Sw(gauss) > p 2. X ■ (54) 



5o(T)r2(m) V^ad, 

For a 100 MW reactor where So(^)^o("^) ~ 1 Tm^, the adiabaticity constraint does not 
impose a serious restriction on the choice of design parameters. 



IV. TRAPPED PARTICLE INSTABILITY 



Now let us consider the influence of the trapped particle mode to the stability of the 
KSTM. The crucial aspect to establish stability of the trapped particle mode is to have 
enough electrons that communicate between the central-cell region and kinetic stabilizer 
region. To determine if stabilization of the fast-growing trapped particle mode is possible, 
we estimate the fraction of electrons in the kinetic stabilizer region that need to reflect back 
into the central plasma [13j. 

To obtain a reasonable estimate requires we need to describe the structure of the potential 
well in the region outside of the plugs. Immediately outside of the MHD destabilizing plug 
(shown on the right-hand side of figure |6| there exists an ambipolar potential energy of 
electrons, "^{B) = — |e|$(i?) (with the choice $vk = 0). This ambipolar potential forms in 
the region between the wall and the plugs, preventing all but the most energetic electrons 
from escaping. The effective potential acting on electrons is Ues{B) = \E'(-B) -|- ^B. Due 



21 



to the opposite sign of their electric charge, the ambipolar potential forces an immediate 
escape of any ions that reach the expander region beyond the plug. 




FIG. 6: Schematic diagram of the negative of the electrostatic potential energy felt by the electrons 
for the focused (solid) and unfocused (dashed) cases. In region Rl, the ambipolar potential is 
determined by the electron temperature of the tandem-mirror-confined plasma. Region R2 is the 
plateau region where the electron density and ambipolar potential are nearly constant in space. 
Region R3 is the kinetic stabilizer region of large stabilizing concave curvature. 

In the region labeled Rl in figure |6| escaping ions experience an increase in kinetic energy, 
E = Tic + i^Tco that is several times the central-cell electron temperature, with k ~ 5. The 
ions in this escaping beam are moving nearly parallel to the magnetic field. We assume that 
there is a rapid transition around a stand-off region where B = Bst- At the stand-off region 
there is a transition in the ambipolar potential; towards the central cell region (to the right 
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of -Bst in figure [6]). The ambipolar potential energy varies on the scale of the central-cell 
electron temperature. Toward the kinetic stabilizer region (to the left of i?st in figure [6]) the 
ambipolar potential tracks with the significantly lower electron temperature of the expander. 
The stand-off position may be determined at the point where the total stress tensor of the 
effluent matches the pressure of the kinetic stabilizer beam. However a proper theory for this 
intuitively described model of the potential structure on either side of the abrupt stand-off 
region still remains for future work. 

The characteristics of region R2 are determined by the kinetic stabilizer beam. As has 
been discussed in section [Tl| in region R2, between the KS region and the stand-off point, 
there exists a density plateau where the ambipolar potential and ion density are nearly 
spatially constant. This plateau region exists regardless of whether the KSTM beam is 
focused or unfocused; however, the plateau region has significantly higher density in the 
unfocused case. 

We assume that the electrons in the kinetic stabilizer have a long mean free path and can 
be described by a Maxwellian distribution with temperature Teks- For simplicity, we take 
Teks <^ ^iEq, although we stretch this inequality to its limits of validity. The Maxwellian 
assumption can be justified for the focused case where an ambipolar potential produces 
a dominant fraction of trapped electrons (shown in figure [6] in region R3 as a solid curve 
between i?w < B < i?ks); the majority of the electrons are trapped and confined long 
enough to relax by collisions to a Maxwellian distribution [3]. For the unfocused case the 
determination of the mean electron energy in the R2 and R3 regions is far more complex. The 
electron distribution may not be a Maxwellian distribution in this case. For the unfocused 
case the ambipolar potential monotonically increases to the wall values. 

For the focused case, the schematic view of the ambipolar potential is shown in the region 
R3 of figure |6| This region is dominated by a large population of trapped electrons that 
neutralize the incoming ion beam but do not readily communicate with electrons in the 
central-cell. The bulk of the electrons in the kinetic stabilizer region are insulated from the 
central-cell plasma by an ambipolar potential. The case with maximum MHD stabilization 
is particularly susceptible to the trapped particle instability [5]. 

The unfocused case leads to moderately reduced MHD stabilization, but has a far less 
restrictive condition for achieving trapped particle mode stabilization. Because the elec- 
trostatic potential monotonically increases in the direction from the expander toward the 
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central-cell, electrons with magnetic moment /i = feel an inwardly-directed electrostatic 
force. However, electrons with a finite magnetic moment feel an outwardly-directed mirror 
force that by itself would produces a magnetic trap what would prevent entrance to the 
central-cell. Nevertheless we presume that if any electron from the distribution of tem- 
perature Teks penetrates through the stand-off point B = Bst, it will also penetrate into 
the central-cell because of the large electrostatic forces directed toward the central-cell for 
B > Bst- For this to happen requires Teks-Bmx/-Bst < i^T^o, where k is a numerical constant 
of order unity. However the spatial potential profile is an important issue that may still 
prevent some electron penetration. For example, if most of the large potential drop is very 
close to the plug, a large fraction of the electrons would be reflected by the magnetic field 
and reduce the trapped electron penetration into the central-cell. Nonetheless, we choose 
the most favorable case for satisfying the trapped particle instability: the case where all 
electrons that reach the abrupt stand-off section penetrate into the central region of the 
plasma. 

Ions moving toward the central-cell from the end wall do not have enough energy to 
penetrate beyond the stand-off region since the large ambipolar potential there repels them. 
These ions are blocked from penetrating by both the repulsive quality of the magnetic 
mirror force (accounted for in the distributions we have chosen) and ultimately by the 
extremely strong electrostatic force in the transition region that confines the hot central-cell 
electrons. In contrast, electrons in the kinetic stabilizer region, whose source is from a supply 
originating from the end-walls, feel forces pushing them into the central-cell. Hence, there 
are no confined ions sharing both the central-cell and kinetic stabilizer region. There are 
electrons that circulate both in the central-cell and in the kinetic stabilizer region. These 
electrons are a source for the charge uncovering term that is a stabilizing factor for the 
trapped particle mode [13]. 

The MHD stability criterion that we use in this paper strictly is the stability criterion for a 
low beta m = 1 flute mode where the radially perturbed displacement is constant everywhere 
along a cylindrical-like column and we shall apply this criterion to the finite beta regime. 
Consider a test function that pinches off the displacement so that there is zero displacement 
in the regions of favorable curvature. The trapped particle mode arises because the pinched- 
off potential removes the stabilizing influence from the good curvature region. In ideal MHD, 
the pinching process excites bending energy that allows the continued stabilization of the 
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mode. However, kinetic processes enable the excitation of an electrostatic perturbation 
that does not excite bending energy. In addition, a low-beta flute mode is electrostatic, 
which enables relatively efficient coupling to other electrostatic perturbations. In a tokamak 
such coupling leads to the Kadomtsev-Pogutse trapped particle mode [H], which allows a 
curvature-driven instability to persist, but at a greatly reduced growth-rate compared with 
the prediction from MHD theory. In a tokamak this reduced growth-rate occurs due to the 
circulation of passing particles through the entire plasma. In a tandem mirror there is a 
substantially lower fraction of connecting particles than in a tokamak (where the fraction 
is near unity). Thus the relative growth-rate of the trapped particle mode is larger relative 
to the trapped particle growth-rate in a tokamak. When the fraction of trapped particles is 
sufficiently low and no additional stabilization mechanism is employed, the growth-rate can 
even be as large as the MHD growth rate of the central-cell and plugs. Mode stabilization 
of the trapped particle mode in a tandem mirror relies on the property that there can be 
a different fraction of connecting electrons to connecting ions. This difference leads to a 
stabilizing effect, which is considered below. 

The stabilization condition for the m = 1 mode does not have the usual finite Larmor 
radius stabilization [15j. However, a similar stabilization mechanism arises due to charge 
uncovering, the difference between the number of electrons and ions that sample both the 
central cell and the kinetic stabilizer regions. In ideal MHD, the cross-field currents due 
to the lowest order E x B drift do not produce an electrical current due to cancellation 
of electron and ion flow velocities. Because of finite Larmor radius effects, the ion E x B 
drift differs from the electron E x B drift. This leads to a current that produces charge 
accumulation. For trapped particle and for electrostatic modes, a similar current emerges 
due to the difference in electron and trapped ion particle populations. 

We analyze the case where density and pressure profiles are Gaussian and have the form 
R{z,iIj) = r{z) exp[—ilj/2iljQ] where ip = iplO) = Borl/2 = B{z)r'^{z)/2. We follow past 
trapped particle mode studies in describing the electrons in the trapping region (in this 
work the KS region) by a Maxwellian distribution at a fixed temperature, Teks- Neglecting 
the effect of rotation found in their work. Berk and Lane[16j found the stability condition 
to be: 

(a;:AQ)^>47^HD(l + Q) , (55) 
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where u* = Teks/lej-Bo'^o- The square of the MHD growth rate 7mhd defined to be 

L 



7mhd = -r • (56) 



L 



dzrQ{z)p^{z) 



where Pra{z) is the mass density. Q is defined 



Q = , (57) 



dzr^{z)pm{z) 



-L 



and 



~2Borl^ J dzrl{z)n^t{z)ej 
oolAQ = . (58) 

dzrQ{z)p^{z) 

-L 



where j denotes the species of particle with charge ej. The density of particles of species j 
at position z that connect is rictjiz), i.e. these particles travel directly between, the kinetic 
stabilizer region and the central-cell when a long mean-free-path electron limit is assumed. 

In the kinetic stabilizer region the eigenmode excitation is assumed to be negligibly small. 
The ambipolar sheath prevents ions outside the central region from returning to the central 
region, so that only electrons connect to the central-cell. Thus in the expression for Q and 
AQ, only electrons contribute to the sums. For this case AQ = Q, so that we find fulfillment 
of the trapped particle stability criterion requires. 



/ Tmhd V 

V7MHD / 



(59) 



If this condition is not met, the growth-rate of the trapped particle mode is approximately 
given by [16] 



7tp ~7mhd/v1 + (5- (60) 
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If Q < 1 the intrinsic growth rate is purely MHD, as if no kinetic stabihzer were present. 



For simphcity, as the stabihty condition for (59), we only use Q > 4(7MHD/i^e) of (59), 
a necessary condition for the stability of the kinetic stabilizer to the trapped particle mode. 
In terms of physically intuitive parameters, the necessary stability condition is 

Te^^ j dzrl{z)n,t{z) > A j dz {p^{z) + p\\{z)) rl{z)'^^-^ . (61) 

ks 

where ks refers to the z integration over a single kinetic stabilizer region. 

We investigate conditions that satisfy both MHD and trapped particle stability. To clarify 



the stability boundary we consider the MHD marginality parameter M^g defined in (39). For 
an MHD stable system 

Lc+L-pi 

dz{p^+p\\) ro — 

= -^^ ^ < 1. (62) 

dz (p± rjj- ° 



dz^ 



ks 



It is convenient to express the trapped particle stability condition relative to the MHD 
stabilization drive of the kinetic stabilizer as: 



Teksy dzrl{z)nct{z) 

ks _ ^eks-^kC"p / '^ct 



4 / dz (p±{z) + pi\{z)) r, 



° dz^ 



> 1 . (63) 







In this expression, {rtct/n) is the fraction of electrons that have orbits in the kinetic stabilizer 
and also penetrate into the central-cell. This fraction is 

dz 



n. 



ct \ ks 



n / I dz . ' ^^^-^ 

ks 
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The numerical factor Ak from the ratio in ( 63 ) is 



dz 



ks 



f dx \ Eq J 



= (65) 

1 



dz I -, --7k(2 — X 



ks 

dz 



5X1/2 



dz- 



B 



(66) 



ks II 



5X1/2 



The last approximate integral in (66) applies when the ion beams satisfy 



A/i5ks/(^o - 5/5t) < 1 . 

For beams with AfiB^/EQ ^ 1, i.e. unfocused beams, we obtain Ak = 1/2. For well-focused 
beams, where St = -Bks and AfiBy^g/Eo ^ 1 we have Ak = 1. 

For a focused beam, the ion density build-up at i? ~ i?x allows for increased MHD 
stabilization of a factor of 1.5-3, as shown in figure |5| However the focusing creates an 
ambipolar potential that prevents most of the electrons within the kinetic stabilizer from 
connecting to the central-cell. For this reason the focused case may have the best MHD 
properties but is extremely susceptible to the trapped particle mode because the factor 
{rict/n), the fraction of trapped and connecting electrons, is exponentially small. 

We therefore examine the unfocused case which will have a substantially larger value of 
{rict/n); for this case it may be possible to satisfy the trapped particle stabilization criterion. 
An unfocused beam results in Ak = 0.5 and the trapped particle stability criterion is 

To apply this stability condition, we have calculated the normalized connecting particle 
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electron density nct{z)ln{z) 



n{z) 



B{z) 



exp 



\B{z) 



1 



(68) 



where iplz) = — |e|0(z)/Toks- The derivation of this formula is given in the Appendix. Recall 
that the subscript st refers to the stand-off position. The connecting fraction of electrons 



{rict/n), defined in (64), is evaluated numerically. The results of this calculation are shown 
in figure [7j 
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FIG. 7: The connecting fraction of electrons in the kinetic stabihzer region that reach the central- 
ceU. The relative radius of the kinetic stabilizer entrance is r\^s/i^w = 0-6. The various curves 
are for different values of e = A/z-SksZ-E-o- The connecting fraction is plotted as a function of the 
stand-off position Xks = ?'ks/^w- 



Despite an increase value of {rict/n) compared to the trapped case, it is still difficult 
to satisfy the stability criterion given by (67) even when Mks is substantially less than 
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unity. This is because all the factors on the left-hand side of (67) are relatively small. For 



example, in figure [7| {rict/n) tends to be small, although it can be close to unity if the 
stand-off position can be made close to the inbound entrance to the kinetic stabilizer region. 
Furthermore, if one attempts to raise the ratio Teks/'^^i-E'o to an order unity quantity, the 
emerging ambipolar potential may prevent ion penetration into kinetic stabilizer region. We 
find in a calculation not discussed here [5] that when a self-consistent attempt to construct 
an ambipolar potential is made, raising Teks/'"ii-Eo ~ 0.3 led to a breakdown of the quasi- 
neutrality condition. Though higher electron temperature solutions should be possible with 
the formation of internal potential jumps on the order of a Debye length, it is likely that 
the maximum ratio Teks/j^j-i-Eo that can give penetration of the ion beam into the kinetic 
stabilizer will remain below unity. 



To satisfy the trapped particle stability criterion given by (67), we need to make the 
MHD stability margin parameter M^s as small as possible; at the same time we must ensure 
that the largest acceptable kinetic stabilizer beam power throughput is less than the alpha 



particle power production of the central-cell. It then follows from (43) that M^s lies in the 
interval 



We fold the trapped particle stability criteria of (67) in with the MHD power constraint 



given in (43). For simultaneous stability to the MHD and trapped particle modes, as well 
compatibility with the power restrictions of the kinetic stabilizer beam, the kinetic stabi- 
lizer's window is 



mm 



/^ct\ T'cks 



< < (70) 



The trapped particle is always more restrictive than the MHD criterion. We take Ac = .175, 
although design improvement from our nominal parameter choice can reduce Ac. Hence the 
compatibility criterion becomes 

^ ^ (Ac/0.175) ^ 35w(g) ^'/V^c(100kev)~ 
4.6 



V3 ro{m)a^ \Bo{T)Eo{keV) J \ A 

m[EQ\ I n 



Bo{T)J V2Tcks; \2.5n, 



< 1 • (71) 
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We have arranged the left hand side of (71), so that each of the bracketed terms have a 
nominal value of unity. We see difficulty in fulfilling this condition. We need experimental 
designs that will allow the bracketed terms to achieve smaller values. Much of the detailed 
physics factors are buried in the factor Ac. Perhaps novel ideas can be developed to gain 
a large reduction factor in this parameter. The nominal value of other factors, such as 
{miEo/2Teks) and (n/2.5nct) have been selected to have as small a value as deemed possible. 
Hence, the present theory for a coUisionless trapped particle mode indicates that there is a 
significant stability/power issue for the stabilization of a symmetric mirror machine with a 
kinetic stabilizer. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the compatibility of the kinetic stabilizer to both MHD stability and 
trapped particle stability. With more beam input power, a smaller value can be obtained 
for the MHD stability parameter, M^s, which must be less than unity to fulfill the MHD 
stability criterion. However, for a fusion device the power requirement for sustaining the 
kinetic stabilizer beam must be substantially less than the fusion power produced. 

In our analysis the nominal maximum beam power is taken to be the fusion alpha power 
production, approximately 20 % of the total fusion power production. Then the MHD and 



power constraints lead to an allowable range of values for M^s, given by (52). In addition 



the window of operation for simultaneous fulfillment of power requirements and trapped 



particle stability leads to the relation given by (71). As it stands, the trapped particle 
instability criterion, together with power constraints, would not be satisfied in a burning 
plasma. Designs that improve on our choices for the of nominal parameters is needed. Below, 
other caveats to this conclusion are discussed. 

A pertinent issue is how severe the trapped particle instability can be. A systematic 
experimental study of this instability has yet to be undertaken. In the kinetic stabilizer 
region, if the electrons are in the short mean-free-path regime, the trapped particle instability 
growth rate is likely to decrease. Then it may be feasible to implement feedback techniques 
to prevent or reduce the harmful effects of the trapped particle instability. Further studies, 
especially experimental studies, are needed to establish a data base to assess the implications 
of exciting the trapped particle instability and determine whether the harmful effects of this 
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instability can be mitigated. 

Our study also indicates the need to develop calculations that can systematically calculate 
the ambipolar potential in the kinetic stabilizer regime. The shape of the electric field in the 
case of unfocused beam propagation is particularly challenging to evaluate, because electrons 
in the kinetic stabilizer are in contact with the wall. In one transit, electrons are unlikely to 
be described by a Maxwellian distribution; it is necessary to verify how reasonable the model 
used in this work is for describing the potential structure in the KS region. A successful 
theory would enable the determination of where the stand-off position would lie. 

There is an additional concern regarding the trapped particle stabilization criterion. If 
trapped particle stability is achieved by improving the communication of electrons from 
the kinetic stabilizer beam to the central cell, there may be a break down in the thermal 
insulation of the hot electrons in the center of the machine. Suppose a value of {rict/n) ~ 0.3 
is achieved, a reasonable estimate of the fraction of connecting electrons necessary to fulfill 
the trapped particle instability criterion. Then the current of cold electron entering the 
central cell will be 30 % of the kinetic stabihzer beam current. These connecting electrons 
are accelerated to the peak energy of the ambipolar potential and replace more energetic 
electrons that leave the system. The energy lost per electron in the exchange is comparable 
to the central-cell electron temperature. The thermal loss rate would be {nct/n)PksTec/ Eq, 
where Tec is the electron temperature in the central cell and Pks is the power sustaining the 
kinetic stabihzer. Since Tec will be very large (> 10^) the connection of the particles, lead 
to an unacceptable power drain. 

Another concern is that the local beta achieved by a focused kinetic stabilizer beam may 
exceed unity. Such a plasma is formed by injecting a beam narrowly distributed in magnetic 
moment so that the beam will reflect back to the wall at designated target field position 
B — Skg. At that position there is large favorable field-line curvature with the magnetic field 
designed so that dr]^/dz — 0. Such a condition leads to a logarithmically large, stabilizing 
MHD response. Mixing the response of a focused beam with a completely unfocused beam 
may enable the best satisfaction of both MHD and trapped particle instability. The increase 
in the MHD response is logarithmic in the small parameter 5 = a/ (ArAyu) / (rks/UT), where 
Tks is the plasma field-line radius aX z = z^s and Ar is the spread in the focusing position 
of the injected particle beam. The local beta increases as Thus if significant MHD 

enhancement is achieved, the local beta value at B^^g is likely to be substantially larger than 
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unity as has been discussed in section HID An issue then arises regarding whether the 
desired focusing can be achieved. 
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Appendix A: Calculation of the density of connecting electrons 

We wish to derive an expression for the density of connecting electrons, nct{z). The 
electrons are assumed to have Maxwellian distribution of the form, 



V Ztc 



(Al) 



where the normalization had been chosen so that the electron kinetic energy and ion potential 
energy ip are in units of Teks/'^^e and that local density is 



oo oo 



n{z) = nw J dv\\ J dv±v±f [v'^) = nwGxp[ijj{z)]. (A2) 







The connecting particles are those that reach the stand-off position z = Zgt with a non- 
zero It follows from energy conservation and magnetic moment conservation that the 
connecting particles satisfy the condition 



7)2 Vu 7)2 R 



Thus 



vl < = ^lmx(^ll) • (A4) 
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The connecting density is then evaluated by integrating over the electron Maxwellian dis- 



tribution inside the kinetic stabilizer with the inequality given by (A4) satisfied. 



nct{z) = nw 



dv 



exp 




dv\\ 
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exp — -+iIj{z) , 



Bst - B{z) 
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Bst - Biz) 

(^st-^(^))i?(^) 
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